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We present the directed flow (vi) measured in Au+Au collisions at % /s JVAr = 62.4 GeV in the 
mid-pseudorapidity region \n\ < 1.3 and in the forward pseudorapidity region 2.5 < \r]\ < 4.0. The 
results are obtained using the three-particle cumulant method, the event plane method with mixed 
harmonics, and for the first time at the Relativistic Heavy Ion Collider (RHIC), the standard method 
with the event plane reconstructed from spectator neutrons. Results from all three methods are in 
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good agreement. Over the pseudorapidity range studied, charged particle directed flow is in the 
direction opposite to that of fragmentation neutrons. 

PACS numbers: 25.75.Ld 



Directed flow in heavy-ion collisions is quantified by 
the first harmonic (vi) in the Fourier expansion of the 
azimuthal distribution of produced particles with respect 
to the reaction plane . It describes collective sideward 
motion of produced particles and nuclear fragments and 
carries information on the very early stages of the colli- 
sion 0] . The shape of v\{y) in the central rapidity region 
is of special interest because it might reveal a signature 
of a possible Quark-Gluon Plasma (QGP) phase O-dQ 

At AGS and SPS energies, V\ versus rapidity is an al- 
most linear function of rapidity Often, just the 
slope of v\ (y) at midrapidity is used to define the strength 
of directed flow. The sign of v\ is by convention defined 
as positive for nucleons in the projectile fragmentation re- 
gion. At these energies, the slope of «i(y) at midrapidity 
is observed to be positive for protons, and significantly 
smaller in magnitude and negative for pions 0, 0, EJ- 
The opposite directed flow of pions is usually explained 
in terms of shadowing by nucleons. At RHIC energies, 
directed flow is predicted to be smaller near midrapidity 
with a weaker dependence on pseudorapidi ty 1 1 CM - . It 
may exhibit a characteristic wiggle |3J, U, LLI1 , whereby 
directed flow changes sign three times outside the beam 
fragmentation regions, in contrast to the observed side- 
ward deflection pattern at lower energies where the sign 
of vi(y) changes only once, at midrapidity. The obser- 
vation of the slope of v\ at midrapidity being negative 
for nucleons or positive for pions would constitute such 
a wiggle. 

In one-fluid dynamical calculations 0, , the wiggle 
structure appears only under the assumption of a QGP 
equation of state, thus becoming a signature of the QGP 
phase transition. Then the wiggle structure is interpreted 
to be a consequence of the expansion of the highly com- 
pressed, disk-shaped system, with the plane of the disk 
initially tilted with respect to the beam direction. Q The 
subsequent system expansion leads to the so-called anti- 
flow |j or third flow component . Such flow can reverse 
the normal pattern of sideward deflection as seen at lower 
energies, and hence can result in either a flatness of V\, 
or a wiggle structure if the expansion is strong enough. 
A similar wiggle structure in nucleon v\ (y) is predicted if 
one assumes strong but incomplete baryon stopping to- 
gether with strong space-momentum correlations caused 
by transverse radial expansion [Jjj- While the predic- 
tions for baryon directed flow are unambiguous in both 
hydrodynamical and transport models, the situation for 
pion directed flow is less clear. RQMD model calcula- 
tions 01 for Au+Au collisions at ^/snn = 200 GeV in- 
dicate that shadowing by protons causes the pions to flow 
mostly with opposite sign to the protons, but somewhat 



diffused due to higher thermal velocities for pions. Simi- 
lar UrQMD calculations predict no wiggle for pions 
in the central rapidity region with a negative slope at 
midrapidity as observed at lower collision energies. 

At RHIC, most of the detectors cover the central rapid- 
ity region where the directed flow signal is small and the 
analysis procedures easily can be confused by azimuthal 
correlations not related to the reaction plane orientation, 
the so-called non-flow effects. Only recently have the 
first u\ results been reported by the STAR Collabora- 
tion |l2j a nd p reliminary results by the PHOBOS Col- 
laboration |l3j . In i the shape of v\ in the region on 
either side of midrapidity is poorly resolved due to large 
statistical errors. This shortcoming arose from having 
only about 70,000 events from the Forward Time Projec- 
tion Chambers (FTPCs) during their commissioning 
in the RHIC run II period (2002). 

In this paper, we present directed flow measurements 
in Au+Au collisions at ^/smn = 62.4 GeV. Results are 
obtained by three different methods, namely, the three- 
particle cumulant method («i{3}) [15j, the eventplane 
method with mixed harmonics (wi{EPi, EP 2 }) [H. Il6j. 
and the standard method Q with the first-order event 
plane reconstructed from neutral fragments of the inci- 
dent beams (ui{ZDC-SMD}). According to the standard 
method 0, directed flow can be evaluated by 

ui{Standard} = (cos(0 - *i))/Res(*i) (1) 

where <p and \&i denote the azimuthal angle of the par- 
ticle and the first-order event plane, respectively, and 
Rcs('I'i) = (cos(4 r i — ^rp)) represents the resolution of 
the first-order event plane. In the three-particle cumu- 
lant method [15|, the flow is evaluated from 



({cos((/) a + <p b 



- 2(f> c ))) = <COS(0„ + <p b - 2(p c )) 
-(C0s((/> a + 0f,))(cos(-20 c )) 

-(cos0 o )(cos(0 6 - 20 c )) 
-(cos0 6 )(cos(0 Q - 20 c )) 

+2(cOS (pa) (COS (pb) (COS(-20 C )) 
= Vl, a Vl,bV2,c 



(2) 



where on the r.h.s. of the first equality, the first term is 
a three-particle correlation and the other terms are to 
isolate the genuine three-particle correlation from spuri- 
ous correlations induced by detector effects. Subscripts 
a, b and c denote three different particles. This method 
was used in the first v\ publication at RHICjl2(. The 
event plane method with mixed harmonics [ijj utilizes 
the second-order event plane from the TPC, ^J PC , and 
two first-order event planes from random subevents in 
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the FTPCs, * FTPCl and * FTPC2 . It measures 
w 1 {EP 1 ,EP 2 } = 



(3) 



^FTPC 



)) 



FTPCi 



FTPC, 



Res(*TPC) 



where the emission angle of the particle <fi is correlated 
with the first-order event plane v]/ FTPC of the random 
subevent (made of tracks of both FTPCs) which does 
not contain the particle. Res(^I / 2 PC ) represents the res- 
olution of the second-order event plane measured in the 
TPC Both the cumulant method and the event 

plane method with mixed harmonics offer enhanced sup- 
pression of non-flow effects, including correlations due 
to momentum conservation, compared with the standard 
method (in which the event plane is reconstructed for the 
same harmonics and in the same rapidity region where 
the event anisotropy is measured). In the present study, 
the procedures to obtain v-j {3 1 and «i{EPi, EP2} are es- 
sentially the same as in Ref. |l2j . In the third method, 



Centrality 


Event plane resolution 


70% 


-80% 


0.179 ±0.005 


60% 


-70% 


0.185 ±0.004 


50% 


-60% 


0.176 ±0.005 


40% 


-50% 


0.167 ±0.005 


30% 


-40% 


0.138 ±0.006 


20% 


-30% 


0.110 ±0.008 


10% 


-20% 


0.081 ±0.010 



TABLE I: The resolution of the first-order event plane Q 
provided by the ZDC-SMDs, as determined from the sub- 
event correlation between east and west SMDs. The errors in 
the table are statistical. 

the reaction plane was determined from the sideward de- 
flection of spectator neutrons (" bounce-off" ) M measured 
in the Zero Degree Calorimeters (ZDCs) |l7(. This is 
the first report from RHIC of flow results with the event 
plane reconstructed from spectator fragments. Five mil- 
lion minimum-bias events were used in this study for each 
of the three analyses, and all the errors presented are sta- 
tistical. Cuts used in the TPC < 1.3) and FTPC 
(2.5 < |7y| < 4.0) analyses are the same as listed in Table 
II of Ref. 0|, except that the vertex z cut is from —30 
to 30 cm. The centrality definition, based on the raw 
charged particle TPC multiplicity with |?7| < 0.5, is the 
same as used previously [Tij . 

In the Fall of 2003, STAR installed Shower Maximum 
Detectors (SMDs) sandwiched between the first and sec- 
ond modules of each existing STAR ZDC at |?7| > 6.3. 
Each SMD consists of two plastic scintillator planes, one 
of 7 vertical slats and another of 8 horizontal slats. The 
two SMDs provide evcnt-by-event information on the 
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FIG. 1: (color online) Directed flow of charged particles as 
a function of pseudorapidity, for centrality 10%-70%. The 
arrows in the upper panel indicate the direction of flow for 
spectator neutrons. The arrow positions on the pseudorapid- 
ity axis corresponds to where the incident ions would lie on a 
rapidity scale. The lower panel shows the mid-pseudorapidity 
region in more detail. The plotted errors are statistical only, 
and systematic effects are discussed in the text. 



transverse distribution of energy deposition of the spec- 
tator neutrons. The weighted center of the energy dis- 
tribution determines the event plane vector on each side. 
The combination of the east and west event plane vectors 
provides the full event plane and the event plane resolu- 
tion is obtained from the correlation of the two event 
plane vectors in the standard way. The «i{ZDC-SMD} 
should have minimal contribution from non-flow effects 
due to the large rapidity gap between the spectator neu- 
trons used to establish the reaction plane and the rapid- 
ity region where the measurements were performed. The 
resolution, as defined in 0] , of the first-order event plane 
reconstructed using the ZDC-SMDs is listed in Table [I] 

The centrality ranges of Au ± Au collisions at ^snn = 
62.4 GeV where the three v\ methods work are slightly 
different: «i{3} fails at centralities less than 5% and cen- 
tralities greater than 70%, because the four particle cu- 
mulant «2{4}, which is a necessary ingredient for mea- 
suring wi{3}, is not measurable in those regions possibly 
due to large vi fluctuations; z;i{ZDC-SMD} fails for cen- 
trality less than 10% because of insufficient event plane 
resolution in central collisions. Figure ^ shows charged 
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FIG. 2: (color online) vi versus rapidity for protons and pi- 
ons. The charged particle v\(rf) is plotted as a reference. The 
different upper end of the pt range for protons and pions is due 
to different limits of the dE/dx identification method. The 
solid and dashed lines are results from linear fits described 
in the text. All results are from analyses using the reaction 
plane reconstructed by the ZDC-SMD, «i{ZDC-SMD}. The 
plotted errors are statistical only, and systematic effects are 
discussed in the text. 



particle v\ as a function of pscudorapidity, rj, for central- 
ity 10%-70% where all three methods work, from Au + 
Au collisions at ^/s^n = 62.4 GeV. The arrows in the 
upper panel indicate the direction of flow for spectator 
neutrons as determined from the ZDC-SMDs. The lower 
panel shows on expanded scales the mid-pscudorapidity 
region measured by the STAR TPC. The results from the 
three different methods agree with each other very well. 
In Ref. the relative systematic uncertainty in vi{3} 
and wi{EPi, EP 2 } was estimated to be about 20%. That 
error estimate was obtained under the assumption that 
the directed flow measurements using two-particle corre- 
lations were totally dominated by non-flow effects. Such 
an assumption provides an upper limit on the systematic 
errors. Ref. [lfj provides further discussion on the sys- 
tematic uncertainties. The comparison of «i{ZDC-SMD} 
and fi{3} indeed shows that the relative difference is no 
more than 20% around mid-pseudorapidity (where the 
directed flow itself is less than 0.005) and the difference 
is only about 5% in the forward pscudorapidity region. 

ZDC-SMD} was also calculated using the information 
from the east and west ZDCs separately as well as sep- 
arately from correlations in the vertical and horizontal 
directions (note that the ZDC-SMDs have a rectangu- 
lar shape); all the results agree within 15%. In another 
systematic study of vi{ZDC-SMD}, a tighter distance of 
the closest approach (dca) cut was applied to reduce the 
number of weak decay tracks or secondary interactions. 
The ratio of v\ obtained with dca < 1 cm to the v\ result 
with the default cut (dca < 2 cm) was measured to be 

dca<lcm / dca<2cm 



1.00 ± 0.07 for charged particles. 



AMPT [3, RQMD 0, and UrQMD [H model cal- 
culations for the same centrality of Au + Au collisions at 
\/ s nn = 62.4 GeV are also shown in Fig. Most trans- 
port models, including AMPT, RQMD and UrQMD, un- 
derpredict elliptic flow (^2) at RHIC energies, and we now 
report that they also underpredict the charged particle 
v\ (77) within a unit or so of mid-pseudorapidity, but then 
come into good agreement with the data over the region 
2.5 < \r}\ < 4.0. While the magnitude of v\ for charged 
particles increases with the magnitude of pscudorapidity 
below |?/| ~ 3.8 for centralities between 10% and 70%, 
our results arc compatible with the peak in \v\\ lying in 
the 1 77 1 region predicted by all three models, namely, 3.5 
to 4.0. 

No apparent wiggle structure, as discussed above, is 
observed within our acceptance. Throughout our pseudo- 
rapidity acceptance, charged particles on a given side of 
7] = flow in the opposite direction to the fragmentation 
neutrons on that side. This is consistent with the direc- 
tion expected in the "anti-flow" scenario but it is also 
the same direction as measured for pions at lower energies 
that is usually related to the pion shadowing by nucleons. 
Assuming that the charged particle flow at beam rapidity 
is dominated by protons, one would conclude that over 
the entire pscudorapidity range v\ (77) changes sign three 
times. However, this does not prove the existence of the 
wiggle structure for protons and pions separately. Mea- 
surements of directed flow of identified particles could 
be more informative in this respect. In STAR, particle 
identification is feasible only in the main TPC, which cov- 
ers the pscudorapidity region |?/| < 1.3. In this region, 
the RQMD model predicts very flat v\ (rj) for pions and a 
clear wiggle structure, with negative slope dv\jdr\ at mid- 
pseudorapidity for protons at ^/s^n = 62.4 GeV. (The 
relatively strong wiggle for pions reported in Ref. jlfj| is 
developed only at higher collision energies.) To maxi- 
mize the magnitude of the possible slope, we select the 
centrality interval 40% to 70%, where flow anisotropics 
normally are at their peak. The result is shown in Fig. [21 
With the current statistics, we observe that pion flow is 
very similar to that of charged particles, with the slope 
at midrapidity dvi/dy about — 0.0074± 0.0010, obtained 
from a linear fit over the region \y\ < 1.3 (dashed line). 
For protons, the slope dv\jdy is —0.025 ± 0.011 from a 
linear fit in \y\ < 0.6 (solid line). At present, STAR's 
statistics for baryons are rather small compared with the 
statistics for all charged particles, and our best estimates 
of the fitted slope are such that a negative baryon slope 
with comparable magnitude to the RQMD prediction is 
not decisively ruled out. For the identified particles, the 
influence of the particle identification procedures on the 
flow values for pions and protons may be a source of er- 
rors. By default we eliminate particles 3<r away from 
the expected TPC energy loss for the relevant particle 
type. When we tightened the cut to 2a instead of 3er, 
we found that for 40%-70% most central events, the 
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FIG. 3: (color online) Directed flow of charged particles as a function of pseudorapidity for different centralities. The plotted 
errors are statistical. 



«i{ZDC-SMD} of pions is reduced by less than 10% while 
the proton ^{ZDC-SMD} stays constant within errors. 

Figure |3 shows vi of charged particles as a function of 
T] for different centralities. We do not observe an onset of 
any special feature in the pseudorapidity dependence of 
v\ at any ccntrality. Preliminary v\ (77) results from PHO- 
BOS [21| for centrality to 40% are consistent with our 
data at the same centrality except that |t>i (77) | from PHO- 
BOS has its peak at \r)\ of about 3 to 3.5, while STAR's 
I vi (77) I peaks at \r/\ about 3.8 or higher. A significant 
change in particle abundances below STAR's transverse 
momentum acceptance cut (0.15 GeV/c) might account 
for some or all of this difference in the \vi \ peak position. 

The transverse-momentum dependence of v\ is shown 
in Fig. 0] Since v\(ri,pt) is asymmetric about 77 = 0, 
the integral of vi(r],p t ) over a symmetric 77 range goes to 
zero. We change v\ (r),pt) of particles with negative 77 into 
— v\(— rj,pt), and integrate over all rj. Due to the small 
magnitude of the i>i signal close to mid-pseudorapidity 
(| 77 1 < 1.3), only the averaged vi(pt) over centralities 
10%-70% is shown. For 2.5 < \r/\ < 4.0, the vi signal 
is large enough to be resolved for different centrality re- 



gions. The poor p t resolution for higher p t in FTPCs lim- 
its the p t range to below 1 GeV/c for 2.5 < \r)\ < 4.0. For 
all centralities and pseudorapidity regions, the magnitude 
of vi is observed to reach its maximum at pt w 1 GeV/c 
for |?7| < 1.3 and at p t w 0.5 GeV/c for 2.5 < |r/| < 4.0. 



Centrality 


Impact parameter (fm) 


70% - 80% 


12.82 + 0.62 - 0.67 


60% - 70% 


11.89 + 0.67- 0.52 


50% - 60% 


10.95 + 0.58 - 0.52 


40% - 50% 


9.91 + 0.47-0.42 


30% - 40% 


8.71 + 0.52-0.31 


20% - 30% 


7.36 + 0.47- 0.26 


10% - 20% 


5.72 + 0.32 - 0.21 


5% - 10% 


4.08 + 0.16-0.21 


0-5% 


2.24 + 0.07- 0.14 



TABLE II: The correspondence between centrality and im- 
pact parameter. 

Note that from its definition, v\(pt) must approach zero 
as pt approaches zero. The centrality dependence of pt- 
integrated v\ is shown in Fig. [5] The values of the impact 
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FIG. 4: (color online) The upper panel shows tji{3} versus pt 
measured in the main TPC (\rj\ < 1.3), for centrality 10%- 
70%. The lower panel shows «i{3} versus pt measured by 
the Forward TPC (2.5 < |r?| < 4.0), for different centralities. 
Note the different scales on both axes for the two panels. The 
differential directed flow of particles with negative 77 has been 
changed in sign as stated in the text. The plotted errors are 
statistical. 



parameter were obtained using a Monte Carlo Glauber 
calculation [2^], listed in Table ITT1 As expected, v\ de- 
creases with centrality. It is seen that v\ in the more for- 
ward pseudorapidity region 2.5 < \rj\ < 4.0 varies more 
strongly with centrality than in the region closer to mid- 
pseudorapidity (\ri\ < 1.3). 

It has been observed that particle emission (both spec- 
tra and flow) as a function of rapidity distance from 
beam rapidity appears unchanged over a wide range of 
beam energies 

El III El, 

a pattern known as limiting 
fragmentation [25J . Figure [5] presents v\ results in the 
projectile frame for three beam energies. In this frame, 
zero on the horizontal axis corresponds to beam rapid- 
ity for each of the three beam energies. The data sup- 
port the limiting fragmentation hypothesis in the region 

-2<y- J/boam < -I- 

In summary, we have presented the first measurements 
of charged particle directed flow in Au+Au collisions at 
v /s - ^7 = 62.4 GeV. The analysis has been performed us- 
ing three different methods and the results agree very 
well with each other. One of the methods involves the 
determination of the reaction plane from the bounce-off 
of fragmentation neutrons, the first measurement of this 
type at RHIC. This method provides measurements of di- 
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FIG. 5: (color online) Directed flow of charged particles as 
a function of impact parameter for the mid-pseudorapidity 
region (\rj\ < 1.3, with the left vertical scale) and the forward 
pseudorapidity region (2.5 < \rj\ < 4.0, with the right vertical 
scale) . The differential directed flow of particles with negative 
77 has been changed in sign as stated in the text. The plotted 
errors are statistical. 
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FIG. 6: (color online) Charged particle vi for Au+Au colli- 
sions (10%-70%) at 200 GeV Gj| (open stars) and 62.4 GeV 
(solid stars), as a function of 77— ybeam- Also shown are results 
from NA49 (circles) for pions from 158^4 GeV midcentral 
(12.5%-33.5%) Pb+Pb collisions as a function of y — 2/beam- 
The 62.4 GeV and 200 GeV points are averaged over the pos- 
itive and negative rapidity regions. All results are from anal- 
yses involving three-particle cumulants, ui{3}. The plotted 
errors are statistical. 



rectcd flow that are expected to have negligible system- 
atic uncertainty arising from non-flow effects. In addi- 
tion, these measurements provide a determination of the 
sign of V\ . In this way, we conclude that charged particles 
in the pseudorapidity region covered by the STAR TPC 
and FTPCs (up to |?y| = 4.0) flow in the opposite direc- 
tion to the fragmentation nucleons with the same sign of 
77. The pi-dcpcndcncc of v\ saturates above pt ~ lGeV/c 
in the mid-pseudorapidity region and pt ~ 0.5GeV/c 
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in the forward pseudo-rapidity region. Over the pseu- 
dorapidity range studied, no sign change in the slope of 
charged-particle v\ versus pseudorapidity is observed at 
any centrality. The centrality dependence of V\ in the 
region of 2.5 < \rj\ < 4.0 is found to be stronger than 
what is observed closer to mid-pseudorapidity. The ra- 
pidity dependence of v± provides further support for the 
limiting fragmentation picture. 
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